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' Coulomb breakup of a projectile in the Coulomb field of a fully stripped heavy nucleus 

' is at present one of the most popular experimental methods to obtain information on 

, reactions of interest in nuclear astrophysics. Its theoretical interpretation presents, how- 

ever, considerable difficulties, due to the three-body nature and the infinite range of 
the Coulomb forces involved. Among the uncertainties affecting present analyses, the 
ki " possible modification of the dissociation cross section by three-body Coulomb final-state 

KN , interactions plays a major role. Various methods which have been proposed to deal with 

■ it are briefly reviewed. However, none of them is based on a consistent and mathemati- 

cally satisfactory quantum mechanical treatment, with the exception of the one proposed 
recently. The latter, being based on the prior form of the dissociation amplitude, makes 
use of a genuine three-charged particle wave function for the final state which is an exact 
solution of the three-body Schrodinger equation, albeit only asymptotically, i.e., for large 
distances. Nevertheless, interesting conclusion can be drawn concerning the influence of 
three-body Coulomb final-state interactions on quantities of astrophysical interest. 
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1. Introduction 

Coulomb breakup reactions, being three- (or more-) problems, are dynamically very 
complicated due to the influence of strong interactions acting between the colliding 
nuclei; and the long-range character of the Coulomb interaction adds non-trivially to 
the complication. However, in the special case of very peripheral collisions the nuclei 
do not come close to each other so that the short-range nuclear forces between them 
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can be neglected. Then the reaction mechanism simplifies considerably: the nuclei 
interact with each other through the electromagnetic field only and unambiguous 
conclusions can be drawn from such experiments. Consider the dissociation 

a + A^b + c + A, (1) 

where particle a is the bound state a — (be), and A is a heavy target. The basic 
mechanism underlying Coulomb breakup is as follows: the fast projectile a moves in 
the Coulomb field of the target A, with the impact parameter being larger that the 
sum of the radii of the nuclei a and A. The electromagnetic target field then lifts the 
fragments b and c into the continuum. There they interact with each other via the 
sum of the Coulomb plus nuclear potentials, but with the target only Coulombically. 

Thus, by choosing appropriate kinematic conditions Coulomb breakup represents 
a powerful, though indirect, tool to obtain information on the radiative capture cross 
section b + c — > a + 7, at astrophysically relevant energies^ (we also recommend the 
reader to consult Ref. El where the general theory of Coulomb excitation and dissoci- 
ation and their applications are reviewed). Its advant ages are exploited around the 
world in different nuclear centers (GSI, KfK, Germany) ^! G ANI L , FranceplNSCL, 
MSU, Cyclotron Institute, Texas A&M University, USaPHHED RIKEN, Japan.EHI) 
It is, however, clearly impossible to cover in a short review all experiments where 
Coulomb breakup reactions have been used to obtain astrophysical information. 

Theoretically, Coulomb breakup reactions are studied, e.g. , by post and prior 
forms of the Distorted- Wave Born Approximation (DWBA), * * ° ' the semi- 
classical approach^] the diffractive dissociati on appro ach and the continuum- 
discretized coupled channel (CDCC) method j ™ * l each of which has specific 
advantages and disadvantages. It must, however, be kept in mind that the astro- 
physical factors ( "5- factors" ) extracted from the analysis of Coulomb breakup data 
are affected by various uncertainties due to (i) the unavoidable extrapolation of the 
data to astrophysically effective energies, (ii) contributions from various electromag- 
netic multipoles (El, E2, Ml), (iii) assum ptions about th e nuclear interactions, and 
(iv) various higher-order effects (see Refs. Hij PU | 21 | 22 |2SI anc [ re f e rences therein). 

One of the most interesting higher-order effects is the so-called post-decay ac- 
celeration (PDA) of the fragments in the Coulomb field of the target. It results 
from the fact that, if the fragments of the projectile have different charge-over-mass 
ratios, the repulsive force exerted by the target nucleus will give rise to a different 
acceleration of each fragment (PDA) which eventually leads to a change in their 
relative energy on their way from the breakup point to the detector. Consequently, 
the experimentally accessible amplitude of the projectile photodissociation will dif- 
fer from the physically interesting one. One of the first direct measurements^ of 
the PDA has been performed for the Coulomb breakup 7 Li — > a + t on 208 Pb, at 
a projectile incident energy E\ = 63 MeV. There it has been been found that the 
minimum of the coincidence spectrum of a-particles and tritons which is expected 
to occur at zero relative energy E at , was shifted due to the different acceleration of 
each fragment in the Coulomb field of the target. Evidently, this different PDA of 
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the a and the triton is due to their different charge-over-mass ratios. The shift of 
the minimum in the coincidence spectrum can serve as an nuclear clock.^1 

Theoretically the PDA effect has been studied by several authors. The sudden 
approximation, which takes into account the interaction to all orders but is valid 
only if the time ordering can be neglect ed, provides a simple method to deter- 
mine the higher-order effects. Calculations^ show that cross sections tend to be 
smaller than those obtained in first- and second-order perturbation theory, at very 
small excitation energy which is the region of interest for nuclear astrophysics. A 
more advanced approach which also takes into account the interaction to all or- 
ders is provide d ^y solving directly the time-dependent Schrodinger equation (Refs. 
| 2U | 22 | 2£i | 2tj | 27 | an( j references therein) . Here, the motion of the projectile and of the 
fragments in the Coulomb field of a heavy target are treated quasi-classically as- 
suming straight-line trajectories, but the relative motion of the fragments quantum- 
mechanically. However, the time evolution in a Coulomb field does not proceed along 
a straight line even for large times because of an additional, logarithmically increas- 
ing time dependence (which is easily seen by i nteg rating the classical equation of 
motion in a Coulomb potential). A recent chedi^of the assumption that the pro- 
jectile and the target do not overlap during the collision ("far-field approximation") 
has shown that it decreases the S'-factor extracted from the data of Ref.Elby 5%. 
However, it is not clear how accurate the quasiclassical approach is in the region 
where the projectile and the target do overla p. 

The CDCC method has also been appliedEd 

As usual the effects of the E2 
transition and higher order excitations are included. Since the relative motion of 
projectile and target, with the huge number of partial waves involved, is treated 
quantum-mcchanically, such kind of approach becomes numerically very involved. 
In some casea^ 1 the semiclassical approximation was seen to compare well with 
corresponding fully quantal calculations at sub-barrier energies. This represents 
a chance for an enormous simplification in the numerical work since it allows to 
check and to possibly supersede the CDCC calculations. At higher energies, where 
also nuclear (grazing) collisions can become important, a Glauber approach is very 
useful!^ Alternat ively , in the stationary framework use has been made of, among 
others, the DWBAl^J and of the quasi-classical approachP2 However, in these in- 
vestigations the three-body scattering wave function for the initial (or final) state 
is approximated by a product of two two-body or quasi-two-body wave functions 
which, as will be explained below, actually completely eliminates PDA. For instance, 
in Refs.^-LG21a "i oca j" momentum is introduced which should allow the incorpora- 
tion of PDA; but in reality it is a "global" momentum because the distance between 
the spectator and the other two particles which are close to each other is rather ar- 
bitrarily fixed at a global value of 10 fm. Besides, the method used provokes further 
questions since the direction of the introduced "local" momentum is undetermined. 

In a systematic approach, the dissociation amplitude should be calculated by 
taking into acc ount the coupling of the different channels, utilizing, for instance, the 
FRESCO code ISE2I However, being presently interested only in an estimation of the 
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possible effects of three-body Coulomb final-state interactions which are responsible 
for PDA on the breakup cross section at small scattering angles and relatively high 
energies, we take resort to some approximations which are admissible under such 
conditions. Starting from the prior form of the dissociation amplitude we introduce 
for the exact three-charged particle wave function occurring therein an approxima- 
tion which, while simplifying the calculations, still preserves the genuine three-body 
nature. Thus, the PDA effect will be accounted for in a systematic manner. 

The plan of the paper is as follows. Section [21 contains a brief discussion why 
Coulomb breakup experiments are useful for obtaining information on capture re- 
actions, and a short summary of pertaining experiments and their astrophysical 
relevance. We then introduce in Sec. 01 the prior-form dissociation amplitude and 
the (approximate) wave function used to describe the three-charged particle final 
state. There we also point out that for an investigation of such breakup processes 
the post form is inappropriate for the derivation of physically correct DWBA or sim- 
ilar approximations. In Sec. 21 some cross section results for the Coulomb breakup 
reactions 6 Li + 208 Pb — ► a + d + 208 Pb and 8 B + 208 Pb — ► 7 Bc+p + 208 Pb are pre- 
sented from which the impact of PDA on the astrophysical 5-factor as determined 
from this reaction can be estimated. Finally, Sec. contains some conclusions. 

We use the units in which h = c = 1. Moreover, unit vectors are denoted by a 
hat, e. g., x = x/x. 

2. Coulomb breakup reactions and their relation to problems of 
nuclear astrophysics 

In nuclear astrophysics, radiative capture reactions of the type 

+ + 7 (2) 

play a very important role. They can also be studied by the time-reversed reaction 

7 + a — ► b + c , (3) 

at least in those cases where the nucleus a is in its ground state. The cross sections 
of reactions J3J and © are related by detailed balance, 

9("9i' 4- \ \ fc 2 

i'^'^' mIw;-) ^ ■ (4) 

The relative momentum between particles b and c is given by k^ c = y/2fn c Ei, c , with 
Hbc being the appropriate reduced mass; the photon momentum is fc 7 = (Eb c + Q), 
where Q is the Q- value of the capture reaction (J2J. Typically one has k^/kb c *C 1 
for energies not in immediate neighborhood of the threshold region kb c ss 0. Due 
to this phase space factor the cross section for the capture reaction (J2J is strongly 
suppressed in comparison to that of the dissociation reaction (J3J . 

In non-resonant charged particle reactions the energy dependence of the cross 
section for reaction is dominated by the penetration of the Coulomb barrier. This 
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energy dependence is usually factored out by denning the astrophysical S-factor 

S(E bc ) = E bc a(b + c -► a + 7) exp(27r??), (5) 

where r\ = Z b Z c e 2 iibc/kbc is the Coulomb parameter. It is the most characteristic 
and important quantity for nuclear astrophysics. 

In fast peripheral collisions of light nuclei with a heavy nucleus the Coulomb 
field of the target nucleus creates the equivalent photon spectrum which allows the 
extraction of the cross section of reaction (J2J) from Coulomb breakup reactionsPlle- 
cent developments, especially of radioactive beams, often permit one to obtain equiv- 
alent information with higher-energy beams. The high-energy Coulomb breakup ex- 
periments generally yield higher event rates; they also yield sometimes information 
which is not available in the direct radiative capture reactions. 

The 6 Li Coulomb dissociation into a + d has been a test case of the Coulomb 
breakup methodPES] The experiments were carried out with a 156 MeV 6 Li beam 
in Karlsruhe^, and with a 60 MeV beam in Heidelberg!^ This reaction is of as- 
trophysical importance since the d(a,7) 6 Li radiative capture is the only process by 
which 6 Li is produced in the standard primordial nucleosynthesis models. 

The radiative capture reaction t(a, 7) 7 Li is astrophysically relevant for the pro- 
duction of the nuclide 7 Li in the early universe and during the nucleosynthesis 
in s tars . The most accurate direct capture measurement has been performed in 
Ref. GUI. The Coulomb dissociation method to investigate the 7 Li — > a + t has been 
applied, too, under various conditions. We mention here the direct breakup of 7 Li 
scattered at 70 MeV off a 120 Sn target,^ a t 63 MeV off 58 Ni, 120 Sn, 144 Sm, and 
208 Pb targetsj^ and at 54 MeV off 197 Au.^' A rather detailed measurement was 
performed recently, with improved experimental techniques!*^ 

The reaction 9 Li — > 8 Li + n was studied with the Coulomb dissociation of a 
9 Li beam of 28.53 A- MeV at MStPl to get information on the cross section for 
the radiative capture reaction 8 Li(n,7) 9 Li, which is of importance for the nucle- 
osynthesis in inhomogeneous big bang models and in Type II supernovaeP^^ It 
is also important for determining the different primordial abundances of Li, Be, 
B, and C, in order to limit the degrees of freedom in inhomogeneous big bang 
nucleosynthesis.^"^ 

The 14 C(n,7) 15 C reaction plays a role in neutron-induced CNO cycles of stel- 
lar evolution phases beyond the main sequence. It is also involved in the chain of 
reactions in primordial nucleosynthesis in the neutron-rich environment of an inho- 
mogeneous big bang. Therefore, at MSU a beam of 15 C ions of E\ = 35 A- MeV was 
used to study the 15 C -> 14 C + n breakup on C, Al, Zn, S n and Pb targets!^ 

The Coulomb breakup of 14 was measured at 

RIKEIs |43144| and GANIlJSl i n 

order to determine the iS-factor of the 13 N(p, 7) 14 capture reaction which is of 
astrophysical relevance in the hot CNO cycle. 

The reactions 11 C(p, 7) 12 N which plays a role in the hot pp-chain, and 
22 Mg(p, 7) 23 A1 which is relevant for the nucleosynthesis of 22 Na in Ne-ric h no- 
vae, were studied at RIKEN by the Coulomb dissociation methodEHEZEHl Their 
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result^^E3 are consistent with the GANIL measurement El The accuracy of the 
value for the radiative width of the 1.19 MeV level in 12 N was improved. 

Finally, the 7 Be(p,7) 8 B radiative capture reaction is important for the solar 
neutrino problem which has been solved recently after 40 years of discussions and 
measurements by many groups. It determines the production rate of 8 B which is 
the source of high-energy electron neutrinos v e . In fact, the astrophysical S- factor 
for this reaction must be known to ±5% in order that this uncertainty no t be 
the domin ant error in predictions of the solar v e flux. The Kamiokandd^ and 
measurements of the neutrino flux from the sun confirmed the standard 
solar model (cf. Ref. ^j). Also the Coulomb dissociation 8 B — > 7 Bc +p has been 
studied at different radioactive beam facilities under different kinematic conditions 
and energies (46.5, 51.9, 83, 90, and 254 A- MeV) at RIKEnP^S MSuP and 

gsi.EH 



3. General Coulomb dissociation amplitude 

Since in the present investigation the main emphasis lies on the elaboration of 
possible effects of three-body Coulomb final-state interactions on the breakup cross 
section at small scattering angles and relatively high energies, we take resort to 
some approximations which are admissible under such conditions. 

As starting point we choose the exact dissociation amplitude for the reaction 
(|T|) in the prior form, 

j\4 prior) = (^lAVl v a $l +) ), AV = V bA + V cA - U. (6) 

Here, VjA = VjA+VjA i s tne interaction potential of particle j(= b, c) with the target 
particle A containing nuclear {Vj\) and Coulombic (V^) parts. The distorted wave 

<1>| in the initial channel is the continuum solution of the two-particle Schrodinger 
equation for the full optical potential U = U c + U N and describes the scattering 
of particles a and A with relative kinetic energy Ei\ ip a is the a — (be) bound state 
wave function with binding energy Ebc, and ^ ' the full final-state wave function 
for three charged particles b, c, and A in the continuum, with the total energy 

E = Ef + Etc = E\ — Ebc- (7) 

The various kinetic energies are defined as Eb c = k^ c /2fibc and E^ = q? l /2[j, a A', k;, c 
is the relative momentum of particles b and c in the exit channel, q^f the relative 
momentum of a and A in the entry, and of the center of mass of (6 + c) and A in 
the exit channel, respectively. The reduced mass of particles i and j is denoted by 
Hij = rmmj/ (rrii + mj). 

We assume that the incident center-of-mass energy E\ of the projectile a is 
high enough and that the scattering angle 8 of the center of mass of the pair (b + 
c) is so small (at given E{) that the impact parameter is much larger than the 
nuclear interaction radius i?N between a and A. Such conditions are typical for 
Coulomb breakup processes. Then all potentials in © can be approximated by their 
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Coulombic parts, AV « AV C = V b c A + V® A - U c . Note that under such conditions, 
<&; is just a two-body Coulomb scattering wave function for the potential U c . 

The dissociation amplitude Jfjjl requires knowledge of the three-charged particle 
wave function ^[ ^ which is a solution of the three-body Schrodinger equation with 
incoming- wave three-body boundary conditions, 

(E — T r — T p — V bc - V bA - VLO^LU ( r ' P) = °' ( 8 ) 
Here, r = r b — r c and p = (mf,r c + m c r;,) /(m(,+m c ) — r^, where Yj is the coordinate 
vector of particle j, while k;, c and qt are the momenta canonically conjugated to r 
and p, respectively. Hence, T T and T p are the kinetic energy operators corresponding 
to the relative motion of b and c, and to the relative motion of the center of mass 
of the pair (b + c) and the target A, respectively. 

Since such a wave function is neither known and nor numerically calculable at 
present, physically justifiable approximations are unavoidable. Such can be found 
by the following consideration. The amplitude JHJ contains the overlap of the initial 
bound state wave function y a (r) for the bound pair (be) and the three-body (6 + 
c + A) final-state continuum wave function. The former limits the range of variation 
of r. But the range of p is not confined and, hence, extends to infinity. Thus, if 
p ^> R' N , where R' N is the larger of the nuclear interaction radii between particles b 
or c and A, it suffices to find solutions of the Schrodinger equation restricted to the 
asymptotic domain £Ia : {r <ti p, R'^ <C p}, i.e., where the distance r between the 
fragments b and c is much smaller than that between their center of mass and A. 

It is easy to see that in the Schrodinger equation © reduces to 

{E T r T p - V bc (r) - U c (p)} *i2 iqf (r, p) = 0(r 2 /p 2 ), (9) 

since there 

V bA (n) + V cA (v c ) ~ V b C A (r b ) + V c c A (r c ) = U c (p) + 0{r 2 /p 2 ), (10) 

with 

V? A {r b ) = Z J Z A e 2 /\p + e J p bc r/m J \, U c {p) = Z a Z A e 2 /p, j = b,c. (11) 

Here, £j = +1(— 1) for j = 6(c), and Zje is the charge of particle j. 

Since the Hamiltonian in Eq. @ separates in the variables r and p, the obvious 
solution is provided by a product of two two-body scattering wave functions which 
we write for later purposes as 

*|- } (r,p) = < c } , qf (r,p) « e**.-**** ^(r)*£>(p). (12) 

Here, V*k ( r ) = e ikbc ' r '0£~' > (r) describes the relative motion of particles b and c, 
interacting via the potential Vb c (r); and \p) = e lqrp ^q f \p) is the scattering 
wave function of the center of mass of (b + c) and A, interacting via the optical 
potential U (p). When inserted in the exact expression ©, this leads to the DWBA 
amplitude conventionally employed in analyses of Coulomb breakup reactions, 

M5»WBA = ( ^(-)^-) |AF C | ^ $ ( + ) ) . (13) 
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We mention that frequently a wave fu nction \l/q f ■* (p) calculated with the full optical 
potential U c (p) + U N (p) is used here! 11 ! 13 ! 

Approximation l|12|) which has dramatic consequences is motivated by two ar- 
guments: (i) the right-hand side of Eq. (|12|l is indeed an exact solution of the 
Schrodinger equation restricted to the domain fl A which is relevant for Coulomb 
breakup processes of interest here, and (ii) its use greatly simplifies the evaluation 
of the dissociation amplitude ©. 

While the second argument obviously applies (see below) , the first one is grossly 
misleading. Namely, it comes as a surprise that the condition that Eq. EH be satisfied 
up to 0(r 2 1 p 2 ) order does not suffice to provide a unique solution)**} To remove 
this non-uniqueness and to select the physically acceptable solution, an additional 
constraint has to be imposed. Such is provided by the requirement that the correct 
solution of the Schrodinger equation © in fl A has to smoothly match with the 
solution of the exact Schrodinger equation JHJ) in the asymptotic region where 
all three particles b, c, and A, are well separated, i.e., r — > oo, p — > oo, since 
^IoH^Ia 7^ 0. The leading term of the asymptotic solution in the asymptotic domain 
Slo is given by the Redmond three-body distorted wavej^ 

yjj(-) (r p) R§ e ik bc-r+iqeP e ~iVbc In (be e ~iribA In ( bA e -ir] aA In ( aA 

« e* k - r +^ ^;\r bc )^\r bA )^;\r cA ), (15) 

where the second for is asymptotically equivalent to the first one. Here, rjij is 
the Coulomb parameter, Qj = fey r\j — ky • ry, and r y (ky ) the relative coordinate 
(relative momentum) between particles i and j. The wave functions "0 k = 
gikij-rij^cc )( r „) describe the scattering of particles i and j, interacting via the 
Coulomb potential Vy (ry) (with incoming wave boundary condition). 

It is apparent, and in fact a simple exercise to check, that expression (|12|) does 
not match in the leading order with the asymptotic solution l|14l) in f2o and, hence, 
is not admissible. That defect notwithstanding, it is conventionally taken as the 
approximation for the final-state three-body wave function ^[ ' . 

Thus, an alternative approximation for ^ ' must be found which, while satisfy- 
ing the above constraints, nevertheless renders the dissociation amplitude calculable. 
A th ree-bod y wave function which satisfies all these conditions has been derived in 
Refs. | 56 | 59 | leading order terrr!^ can be written as 

* ^• r+iqf - p t ) ( P )( r )C]( r ^)4:](rcA). (16) 

The new wave function ^7 (p) ( r ) ^ e lki>c( - p ) Tba '4 > ]^ > ( r ) describes the scattering 
of particles b and c, interacting with the full potential V{, c (r) and moving with 
relative momentum k;, c (p). As a consequence of the long-range Coulombic three- 
body correlations between b and A and between c and A the latter, however, depends 
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on the distance of their center of mass to the charged particle A, and is given by 

VbA (k b A + P) VcA (k cy i + p) 



k bc (p) = k hc - a(p) := Hbc 



(17) 



m b (l + p.k bA ) m c {l + p-k cA )_ 

It is not difficult to see that the wave function l|16f) has the desired properties: 

(i) it is asymptotically solution of the asymptotic Schrodinger equation @ and 

(ii) it smoothly matches with the Redmond wave function in Qq and, hence, is 
asymptotically also solution of the original Schrodinger equation JSJ. 

It is apparent that the three-charged particle wave function l|16(l which, as we 
emphasize once more, is the only physically admissible, asymptotically exact solu- 
tion of the three-body Schrodinger equation @ for use in the dissociation amplitude 
©, is much more complicated than (the asymptotic part of) the simple factorized 
form (|12f> conventionally used in DWBA calculations. Insertion into © yields for 
the appropriate dissociation amplitude 

M ; ( f prior) » (^[- )(as) |Al/ c | <p a *[ +) ). (18) 

Since the integration over r is limited by the bound state wave function ip a (r) of 
the projectile a, and if the impact parameter is larger than the nuclear interaction 
radius i?N between a and A, we can utilize the multipole expansion for AV. It leads 
to the multipole expansion of the dissociation amplitude l|18fl : 

M « = E wti Mlm - (19) 

LM 

The multipole transition matrix elements M^m are defined as 
M LM = Z A Z L cS) e 2 J dp J dr^- ){as) *(r,p) 

x ~L~+i Y lm (r) Ylm (p) ifa (r) $ W (p) , (20) 

with Z L cS) = [ifa [(-l) L Z b / m L + Z c /m^] and $| +) (p) = Note that if in 

v j,(-)( a& ) a j g0 ^.jjg ne xt-to- leading-order terms are taken into account, two more terms 
of similar structure but suppressed by an additional factor 1/p occur in H20JI . J 9 

The right-hand side of Eq. (|20|l contains a six-dimensional quadrature. This is to 
be contrasted with the conventional approach. There, the multipole matrix elements 
of the DWBA amplitude (fTl^l separate into a product of two terms, each requiring 
three-dimensional quadratures. This greatly reduces the numerical expense. 

Use of the amplitudes (jl8|) . resp. (|2C)|> . which are based on the prior- form dis- 
sociation amplitude has two crucial consequences. First, since the wave function 
occurring therein has the correct asymptotic behavior this expression is consistent 
with the general principles of scattering theory of charged particles. This is to be 
contrasted with the standard DWBA of the dissociation amplitude (| 1 31) the deriva- 
tion of which usually starts from the post form, A/ ; ( f post) = ^aA \ AV \ 
with AV as in Eq. 10. The final-state wave function is given by the right-hand side 



February 5, 2008 14:32 WSPC/INSTRUCTION FILE Coul'sing'subm 



10 



of Eq. (|12(l . The exact initial-state three-body wave function is then approximated 
by the product ^fi + ' w tp a $>{ + \ However, in the presence of Coulomb potentials 
acting between all three particles, the post-form amplitude is not valid, in contrast 
to the prior form and hence is not equivalent to the latter. The reason is that 
the final-state product wave function occurring therein does not produce the correct 
behavior in any asymptotic region of configuration space. Secondly, as mentioned 
earlier, the PDA effect occurs for Zb/m,b ^ Z c /m c , since it is caused by the differ- 
ence in the Coulomb forces exerted by the target on each fragment b and c, which 
leads to the different accelerations. This fact gives rise to a change of the relative 
velocity of the fragments when they finally have reached the detectors as compared 
to their relative velocity at the moment of the decay of the projectile a. Use of 
approximation 1|12() . however, completely eliminates the PDA effect because (i) the 
wave function part describing the scattering of the fragments b and c in the Coulomb 
field of the target A is approximated by \p) which describes the scattering of 
the center of mass of the pair (b + c) off the target rather than the scattering of 
the individual fragments, and, related to it, (ii) the relative motion wave function 
' (r) of the fragments b and c in the continuum contains the asymptotic relative 
momentum between b and c and, hence, does not allow for its modification due to 
the action of the target Coulomb field between the breakup point and the counter. 



4. Applications 

4.1. 6 Li breakup into ot-particle plus deuteron 

As a first application of the proposed method we have computed^ the triple- 
differential cross section (TDCS) for the reaction 

Pb — >a + d + 20S Pb, (21) 



3 Li 



208 



for a projectile energy Ei = 156 MeV and scattering angles 9 a = 6 a = 3°. In 
this investigation only the leading term of the asymptotic wave function (|16fl has 
been retained. Moreover, because of the peripherality of the reaction it suffices 
to take into account only the Coulomb interaction between the fragments in the 
final state so that ip^ \(p)^ 1S smT Ply a Coulomb scattering wave function, albeit 
with local momentum h ac i(p). Since for p > i?N> the additional term &(p)/p in the 
local momentum turned out to be much smaller in magnitude than the asymptotic 
momentum k QC ;, this wave function was expanded in powers of \a(p)\/pk a d- Keeping 
only the leading terms we arrive at (tp(z) the digamma function) 

A{p,k 



A(p, k ad ) 



1 



a(p) 



P 

kad Vad 



iip(l - iriad) 



(22) 



(23) 



Use of l|2*2*)l greatly simplifies the calculation of the amplitudes Mlm since each of 
them separates into a product of two quantities contaning integration over either r 



I 



I 
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or p only. Thus, in this approximation the whole effect of PDA resides in the factor 
A(p, k QC (); or explicitly, for A(p,'k a d) = PDA does not occur. 

Both El and E2 transitions were included. In the isospin formalism, for the 
6 Li — > a + d breakup the El transition is forbidden (AT = 0); thus the main 
contribution comes from E2. Nevertheless, since the isospin rule is violated and the 
charge-over-mass ratio Z a /m a of the a-particle differs only little from that of the 
deuteron, El contributes comparably to E2 for E ac i < 100 keV (see also Ref. 
where the contribution of El to the reaction a + d — > 6 Li + 7 is discussed). 

In Fig.^we present the calculated TDCS. Inspection shows that for \E a J\ < 300 
keV the agreement between the experimental and theoretical results is very good. 
Hereby, the notation is such that positive (negative) values of E a d are associated 
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_1 I 1 1 1 1 1 1 

-0.6 -0.4 -0.2 0.2 0.4 0.6 
Ead(MeV) 

Fig. 1. The TDCS for the reaction 6 Li + 208 Pb > a + d + 208 Pb at a projectile energy Ej = 156 

MeV and scattering angle 8 = 3° of the fragments, as a function of the a — d relative kinetic 
energy i? Q( j. Negative (positive) values correspond to a — d relative velocities v ac i < {v a d > 0). 
Experimental data are from Ref. El 

with "forward" and "backward" emission of the a-particle in the center of mass 
system of the decaying 6 Li nucleus. The astrophysical S'24-factor for E a d = ex- 
tracted therefrom (by linear extrapolation into the energy range E a d < 20 keV) is 
S2i{0) = 1.2 nb-MeV. This value agrees with that obtained^ from a calculation of 
the cross section for the radiative capture d(a,7) 6 Li at ultra-low energies En^, but 
disagrees with the old extrapolated value S^ p (0) = 9.1±1.8 nb-MeV in Ref.^™ Note 
that an earlier analysiJ^of the same data gave a value of 624(0) = 1.7 nb-MeV. 

Fig. [21 illustrates the role of the factor A(p, k Q d), which comprises within 
our approximations all three-body Coulomb final-state interaction effects to the 
cross section, in the energy region \E a d\ > 20 keV (because of our assumption 
\A(p, k ac j)|/p <C 1 the curves cannot be extrapolated to the point E a d — 0). Nev- 
ertheless it is apparent that the ratio of the cross sections calculated with the con- 
ventional asymptotic wave function [A(p, k ac j) = 0] to the one calculated with the 
asymptotically correct wave function, which takes into account three-body Coulomb 
effects [A(p, h a d) 7^ 0], behaves differently in the regions v a d < and v a d > 0. That 
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1.2 
1.1 
1.0 

R 

0.9 
0.8 

-0.6 -0.4 -0.2 0.2 0.4 0.6 
E„ d (MeV) 

Fig. 2. Ratio d l a [ A Znl^ 6Ll/ , d ^°' d/ , d Z ad of the TDCS calculated without (A(p,k ad ) = 0) and 
with (A(p, k ac j) ^ 0) allowance for three-body Coulomb effects, as function of a — d relative kinetic 
energy E ad . Negative and positive values of E ad have the same meaning as in Fig. HI 

is, three-body Coulomb effects, together with the interference of El and E2 tran- 
sitions, give rise to an asymmetry of the cross section ratio with respect to the 
point v ac i = 0. Hence, it is compulsory to accurately take them into account when 
deducing the S'24-factor from the TDCS. 

4.2. 8 B breakup into 7 Be plus proton 

As a further example we have investigate d^^ the breakup reaction 

8 B + 208 Pb^ 7 Be +J3 + 208 Pb, (24) 

at 8 B beam energies E { = 46.5 A- Me\EH and E { = 83 A- MevPtaking into account 
dipole and quadrupole contributions. 7 Be-p relative kinetic energies are varied be- 
tween 100 keV and 1 MeV. Recall that this reaction allows the extraction of the 
astrophysical S'-factor for the important radiative capture reaction 7 Be+p — > 8 B+7. 
The final 7 Be nucleus is assumed to be in the ground state. Note that in the kine- 
matic region of interest, both the binding energy £7 Bep = 0.137 MeV of the proton 
to the 7 Be nucleus as well as the kinetic energy of the relative motion of the proton 
and Be in the final state are much smaller than the collision energy E[. 

When u sing the asymptotic three-body wave function including next-to-leading- 
order terma^ for the final state, the breakup amplitude is given as a sum of three 
terms. The first one arises from the leading term l|18|l (of O(p )) and is expected 
to give the largest contribution. Indeed, it turned out that the cross section is 
completely governed by it. The other terms are suppressed by an additional factor 
1/p in the wave function, contributing less than 1% to the cross section. 

The dipole transition to the Be-p s-wave continuum state is found to dominate 
the amplitude; the contribution of the dipole transition to the c?-wave state is small 
at small relative kinetic energies E7 Bep but increases up to 30% at 1 MeV. Moreover, 
in our approach the quadrupole amplitude depends only on the projectile energy 
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Ei and on the 7 Be-p relative energy -EV Bep but not on the scattering angle 8. Its 
contribution relative to that coming from the dipole amplitude changes from a few 
per cent at small E7 Bep and small 9 to up to 30% near the grazing angle. We mention 
that the same behavior of the quadrupole term was observed in a calculation using 
first- and second-order perturbation theory in the semiclassical approach!^ The 
consequence is that for radiative capture reactions at energies of interest in nuc lear 
astrophysics, the quadrupole amplitude is negligible. For details we refer to Ref . ^1 
Typical results are presented in Tables ^ and [21 Table ^ shows the ratio of the 
double differential cross section (DDCS), calculated with the (asymptotic) three- 
body wave function l|16|) . as compared to the one using the conventional ansatz 
(|12fl for the final state, at incident energies E t = 46.5 A- MeV. The former takes 

Table 1. Ratio R(6) of the DDCS for the reaction 208 Pb( 8 B, 7 Bep) 208 Pb , calculated 
with the asymptotic three-body Coulomb wave function 1161 . to the DDCS calculated 
using the conventional wave function 1121 . as a function of the scattering angle 9, for 
various values of the 7 Be-p relative energy. The projectile energy is Ei = 46.5 A- MeV. 















R(6) 












[MeV] 


1.0° 


1.5° 


2.0° 


2.5° 


3.0° 


3.5° 


4.0° 


4.5° 


5.0° 


5.5° 


6.0° 


0.1 


1.16 


1.26 


1.37 


1.49 


1.63 


1.77 


1.91 


2.04 


2.17 


2.30 


2.41 


0.15 


1.02 


1.04 


1.07 


1.11 


1.14 


1.16 


1.19 


1.21 


1.23 


1.24 


1.25 


0.2 


1.01 


1.01 


1.02 


1.03 


1.03 


1.03 


1.03 


1.03 


1.02 


1.02 


1.01 


0.4 


1.00 


1.00 


0.99 


0.99 


0.98 


0.98 


0.97 


0.96 


0.95 


0.94 


0.92 


0.6 


1.00 


1.00 


1.00 


0.99 


0.99 


0.99 


0.98 


0.98 


0.98 


0.97 


0.97 


0.8 


1.00 


1.00 


1.00 


1.00 


1.00 


1.00 


0.99 


0.99 


0.99 


0.99 


0.99 


1.0 


1.00 


1.00 


1.00 


1.00 


1.00 


1.00 


1.00 


1.00 


1.00 


1.01 


1.01 



into account three-body Coulomb final state interactions and, hence, PDA while 
the latter does not. Thus, deviation of this ratio from unity is a direct measure 
of the influence of PDA. We only mention that a calculation at 83 A- MeV shows 
that PDA effects decrease with increasing projectile energy and relative energy of 
the outgoing fragments. But at small relative energies and large scattering angles, 
PDA is important, the deviation from unity reaching more than 10%. Moreover, 
from Table El we conclude that PDA affects the single differential cross section 
(SDCS) less than the DDCS. This is understandable because PDA is large near the 
grazing angle where the cross section is small; and the integration over the scattering 
angle of the DDCS levels out its contribution. We mention that our approach may 
produce a large error near the grazing angle where the nuclear interaction between 



Table 2. The ratio as in TableQ]but for the SDCS da/dE7 Bep , at £ ; = 46.5 
and 83 A- MeV, obtained by integrating the DDCS over the scattering angle 
9 G [2° - 6°] and [0.1 - 1.77°], respectively. 



E 7Bcp [MeV] 


0.1 


0.15 


0.2 


0.4 


0.6 


0.8 


1.0 


Ratio (Ei = 46.5 A- MeV) 
Ratio (Ei = 83 A- MeV) 


1.90 
1.13 


1.18 
1.02 


1.03 
1.00 


0.97 
1.00 


0.98 
1.00 


1.00 
1.00 


1.00 
1.00 
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the projectile and the target becomes important which has been omitted here. 
From these results the following conclusions can be drawn. 

(i) At 46.5 A- MeV incident energy, the DDCS is enhanced by three-body Coulomb 
final-state interactions only at the lowest relative energies Et-q cp while in the medium 
energy range it is reduced; at the highest energies considered their effect becomes 
negligible. That is, the three-body effects increase the probability of dissociation of 
8 B up to energies of Etq op ~ 200 keV, but start suppressing it at ~ 400 keV. 

(ii) Three-body Coulomb effects become more important as 9 increases. The rather 
small influence at very small 9 is due to the fact that there large impact parameters 
dominate for which the three-body correlations are weak. Thus, in order to avoid 
PDA effects, experiments at E{ = 46.5 A- MeV should be performed at 9 < 4°. 

(iii) On the other hand, PDA effects decrease as the incident energy increases, cf. 
Table which co mp ares the SDCS ratios with and without three-body effects, at 
Ei = 46.5 A- Me\^21and 83 A- MevP Thus, PDA effects can also be made negligi- 
bly small by choosing sufficiently large projectile energies. 

(iv) Since our present approach is based on an asymptotic three-particle wave func- 
tion and, hence, does not allow to perform the calculations for small enough 7 Be-p 
relative kinetic energies, extrapolation to E7 Bep = seems too arguable. Hence, no 
reliable quantitative predictions about the influence of three-body Coulomb effects 
on the Coulomb dissociation cross section and the 5i7-factor to be extracted from it 
can be made. However, qualitative conclusions, in particular under which kinematic 
conditions PDA effects are negligible, can already be drawn. This shortcoming is 
shared also by the standard DWBA approaches based on Eq. I|13|) since there also 
only a solution, and indeed an unacceptable one, of the asymptotic Schrodinger 
equation and not a genuinely calculated three-body wave function is used. 

5. Summary and outlook 

Peripheral collisions of medium- and high-energy nuclei, that is collisions in which 
these nuclei pass a heavy, fully stripped target nucleus at distances beyond the 
nuclear interaction radius and which are therefore dominated by electromagnetic 
interactions, have evolved into an important tool of the nuclear astrophysical re- 
search, because of the intense source of quasi-real photons available from the target 
in such collisions. In addition new experimental facilities, in particular the present 
and forthcoming radioactive beam factories, render feasible the efficient investiga- 
tion of photo-interactions with (stable and unstable) nuclei, such as single- and 
multi-photon excitations and electromagnetic dissociation. 

In Coulomb breakup processes, among the various so-called higher-order effects 
the post-decay acceleration (PDA) of the fragments in the target Coulomb field 
which is due to three-body Coulomb final-state interactions plays an important 
role, especially at low and medium collision energies. Various theoretical and com- 
putational methods presently in use for taking them into account have been briefly 
sketched. We have, in particular, emphasized that the frequently used conventional 
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separation of the dissociation amplitude into a product of a purely two-body exci- 
tation amplitude and a quasi-two-body scattering part of the center of mass of the 
projectile has the consequence that all PDA effects are eliminated. 

To overcome this problem we have proposed an approach which makes use of a 
genuine three-charged particle wave function to describe the final state. The latter is 
solution of the three-body Schrodinger equation, albeit only in an asymptotic region 
of configuration space which, however, is relevant for breakup processes. Because 
of this our estimation of three-body final-state interaction effects is reliable only 
for small scattering angles. For larger angles a more accurate three-body Coulomb 
scattering wave function has to be used which is not yet available. Therefore, further 
investigations of the PDA effects are called for. Nevertheless, we have shown that al- 
ready the use of even an only asymptotically correct wave function for three charged 
particles in the continuum is beneficial for extracting reliable values of astrophys- 
ical factors from Coulomb breakup reactions, or at least can help pinpointing at 
kinematic conditions in which PDA effects should play a negligible role. We finally 
mention that the discovery of a low lying i51-strength in neutron-rich nuclei and the 
demand to determine astrophysical S-factors of radiative capture processes, as for 
7 Be(p, 7) 8 B, inspires new calculations, possibly with new computational techniques, 
of Coulomb breakup processes. Thus, Coulomb dissociation provides comprehensive 
opportunities to study the interaction of exotic nuclei with photons. 

We have not mentioned the Coulomb breakup through excitation of a resonance 
state of the projectile which is expected to lead to interesting results in experimental 
as well as theoretical respect. 

Among the problems still to be solved are (i) from the experimental side the 
accurate determination of the relative energy of the fragments, (ii) from the theo- 
retical side the extrapolation of the astrophysical factors to zero relative energy, the 
accurate determination of the E2 contribution, relativistic corrections and the nu- 
clear interference terms in the cross section which become more and more important 
as the scattering angle increases, especially if it approaches the grazing angle. 
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